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1.  Introduction 

In  many  applications,  nonlinear  equilibrium  problems  typically  involve  a 
number  of  intrinsic  parameters.  Hence  the  resulting  equations  have  the  generic 
form 

(1-1)  H(y,t)  =  y^, 

where  y  and  t  vary  in  some  state  space  Y  and  parameter  space  T,  respectively, 
and  H  is  a  given  mapping  with  domain  in  Y  x  T  and  range  in  Y.  In  general, 
the  set  of  all  solutions  (y,t)  of  (1.1)  forms  a  manifold  in  Y  x  T  and  interest 
centers  on  analyzing  the  characteristic  features  of  this  equilibrium  manifold. 

Usually,  the  parameter  space  T  has  finite  dimension,  btit  the  state  space  Y 
is  infinite-dimensional.  Thus  for  a  computational  analysis,  finite-dimensional 
approximations  of  (1.1)  have  to  be  Introduced,  and  then  methods  are  required  for 
determining  the  features  of  the  solution  manifold  of  the  approximating  finite- 
dimensional  system  of  parametrized  equations.  At  this  time,  the  principal  component 
of  such  methods  always  is  a  general  form  of  continuation  process  for  the  trace  of 
paths  on  the  manifold.  In  addition,  special  procedures  are  available  for  detecting 
and  determining  specific  features  along  these  paths.  The  literature  in  this  area 
is  extensive,  but  there  is  no  need  to  give  specific  citations  since  no  use  of  them 
will  be  made.  It  should  be  noted,  however,  that  there  are  only  few  studies  about 
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the  Integrated  use  of  these  various  techniques  for  a  thorough  analysis  of  equi¬ 
librium  manifolds. 

The  need  for  Introducing  flnlte-dlmenslonal  approximations  of  the  equations 
(1.1)  leads  to  further  questions  about  the  relation  between  the  features  of  the 
solution  manifolds  of  the  original  and  the  discretized  equations  and  the  errors 
Introduced  by  the  approximation.  Relatively  little  has  been  done  so  far  In  this 
area,  and.  In  fact,  the  series  of  articles  [3],  [5j  appears  to  be  the  only  one 
which  addresses  the  a  priori  estimation  of  the  errors  In  the  case  of  a  one- 
dlmenslonal  parameter  space,  and  In  [l]  some  a  posteriori  estimates  of  these  errors 
are  Introduced  for  certain  boundary-value  problems  In  one  space  dimension. 

In  this  paper,  we  outline  some  techniques  for  ensuring  the  existence  of  so¬ 
lution  paths  of  a  rather  general  class  of  nonlinear  equations,  as  well  as  of  their 
flnlte-dlmenslonal  approximations,  and  for  assessing  the  error  between  these  paths. 

A  tool  In  this  analysis  Is  the  theory  of  nonlinear  Fredholm  operators.  The  dis¬ 
cussion  of  the  operator  discretizations  Involve  a  specific  uniformity  condition 
which  represents  a  restriction  of  the  class  of  operators  used.  The  implication  of 
this  condition  Is  the  subject  of  ongoing  research.  But  we  show  here  that  the  results 
apply  directly  to  mildly  nonlinear  elliptic  boundary  value  problems  of  the  type 
considered  in  [3],  [s].  At  the  same  time,  it  appears  that  they  cover  much  more 
general  classes  of  operators  as  well. 
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with  a  fixed  e  F(W).  For  this  purpose,  suppose  that  for  some  1  the 

m-regularity  set 

(2.3)  R  (F)  =  {x  e  W:  dim  ker  DF(x)  =  m,  rge  DF(x)  »  Y} 

m 

of  F  is  nonempty.  Then  the  restriction  of  F  to  R  (F)  is  a  nonlinear  Fredholm 

m 

operator  of  index  m  (see  eg.  [2]).  Our  analysis  of  (2.1)  centers  on  this  fact. 

In  [4]  it  was  shown  that  in  the  case  m  =  1  the  set  Rj^(F)  is  open  if  it  is 
not  empty.  A  slight  modification  of  the  proof  given  there  shows  that  if  ^  'P 

for  some  m  ^  1  then  again  the  set  is  open. 

For  any  x^  e  R(F)  there  exist  closed  subspaces  V  C  X  such  that 
X  =  V  0  ker  DF(x^) .  For  any  such  choice  of  V  the  restriction  DF(x^)|v  is  an 
isomorphism  between  V  and  Y  and  the  inverse 

(2.4)  =  (DF(x^)|v)‘^  e  L(Y,V) 

is  a  right  inverse  of  DF(x^)  on  Y.  With  this  the  equation  (2.2)  may  be  written 
in  the  form 

(2.5)  F(t  +  A^y)  =  Yg.  7  ^  t  e  ker  DF(x^),  t  +  A^y  e  W, 

which  corresponds  to  (1.1)  with  Y  as  the  state  space  and  =  ker  DF(x^)  as  the 
m-dimenslonal  parameter  space. 

The  basic  result  about  the  solution  set  of  (2.2)  may  now  be  phrased  as  follows: 

Theorem  2.1:  Suppose  that  (2.1)  is  given  and  that  nonempty  for  some 

m  >  1.  Then  for  any  fixed  y  e  F(W)  the  regular  solution  set  F^  ^(y  )  0  R  (F) 

—  o  o  m 

of  (2.2)  is  a  nonempty,  open,  m-dlmenslonal  c’^-manlfold  in  X. 

The  proof  is  standard  (see  eg.  a  finite-dimensional  analog  given  in  [8], 
p.  11). 

As  mentioned  in  the  Introduction,  the  computational  procedures  for  analyzing 
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the  solution  manifold  of  (2.2)  consist  principally  of  methods  for  approximating 
paths  on  that  manifold.  Such  a  path  is  defined  as  the  solution  manifold  of  some 
reduced  equation  defined  by  a  Fredholm  operator  of  index  1  on  its  regularity  set. 
Thus  without  restriction  of  the  generality  we  may  assume  that  (2.1)  is  given  and 
Rj^(F)  is  not  empty.  For  ease  of  notation,  we  shall  write  M(y^)  ••  F^  ^\y^)  f\  Rj^(F) 
for  the  open,  one-dimensional  c'^-manlfold  constituting  the  regular  solution  set 
of  (2.2). 

We  consider  first  the  question  of  the  choice  of  suitable  local  parametrlzatlons 
of  M(y^).  It  turns  out  that  a  possible  choice  corresponds  to  a  typical  approach 
used  in  continuation  procedures. 

M(y^) 

ker  DF(x  ) 
o 

Figure  1 

X  +V 
o 

In  line  with  the  earlier  discussion  leading  to  (2.5),  a  local  parametrlzation 
of  W(y^)  at  a  given  point  x^  e  M(y^)  is  defined  as  a  triple  {V,A,z^}  con¬ 
sisting  of  a  closed  subspace  V  of  X,  a  linear  map  A  e  L(Y ,V) ,  and  a  point 

z  e  X  such  that 
o 

(1)  ker  DF(x  )  n  V  =  {0>; 
o 

(11)  A  is  an  isomorphism  from  Y  onto  V; 

(2.6) 

(iii)  z^  ^  V; 

(Iv)  X  ■  V  ©  T  T  ■  span  {z  }. 

O’  o  o 

As  indicated  in  Figure  1,  we  consider  the  family  of  parallel  linear  manifolds 

x  +  tz  +  V,  t  e  IR^. 
o  o 


(2.7) 
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The  conditions  (2.6)  ensure  that  these  manifolds  are  transversal  to  ker  DF(x  ) 

o 

and  hence  that  locally  near  t  =  0  the  solution  manifold  may  be  parametrized 
in  terms  of  t.  This  is  the  content  of  the  following  result: 

Theorem  2.2:  Suppose  that  (2.1)  is  given  and  Rj^(F)  iji.  Moreover,  let  {v,A,z^} 
be  a  local  parametrization  of  M(y^)  at  the  point  e  there  exists 

an  open  interval  J  0  e  J,  an  open  neighborhood  U  c  x  of  x^,  and  a 

unique  C^-map  n:  J  "*■  Y  such  that 

(2.8)  M(y  )niI={xeX:x  =  x  +tz  +  An(t),  t  c  j}. 

O  0  0 

The  proof  follows  from  an  application  of  the  inverse  function  theorem  and 
the  use  of  the  conditions  (2.6). 

In  the  standard  continuation  method,  z^  defines  the  predictor-line  and  the 
corrector  produces  the  step  An(t)  from  the  predicted  point  x^  +  tz^  to  the 
point  x^  +  tz^  +  An(t)  on  M(y^). 

At  any  point  x  e  M(y  )  choose  a  nonzero  vector  z  e  ker  DF(x  ).  Then, 
o  o  o  o 

as  noted  earlier,  there  exist  closed  subspaces  V  c  X  such  that  X  >  V  0  ker  DF(x^) . 

With  any  such  V  and  the  associated  mapping  (2.4),  the  triple  {V,A,.,z  }  is  a 

V  o 

local  parametrization  of  M(y  )  at  x  .  In  the  setting  of  continuation  methods, 

o  o 

these  parametrlzatlons  correspond,  in  essence,  to  the  pseudo-arclength  parametrizatlons 
(see  eg.  [6]). 

The  question  now  arises  how  far  a  local  parametrization  may  be  extended.  For 
this,  note  that  the  set 

(2.9)  A  -  {x  e  U:  DF(x)A  is  an  isomorphism  of  Y  onto  itself} 

is  certainly  nonempty  and  open.  With  this  a  generalization  of  a  result  in  [7] 
may  be  phrased  as  follows: 

1 

1 

Theorem  2,3:  Under  the  conditions  of  Theorem  2.2,  let  C  X  be  the  maximal  con-  j 
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nected  subset  of  M(y  )  H  A  which  contains  x  . 

o  o 

X  r 

terval  J  c  IR  ,  0  e  J  ,  and  a  C  -map  n  :  J  ■+•  Y 

o  o  o  o 


(2.10)  M  =  {x  e  X:  X  =  X  +  tz  +  An  (t) , 

o  o  o  o 


Then  there  exists  an  open  in- 
such  that 

t  e  J  } . 
o 


The  condition  that,  for  x  e  M  ,  the  map  DF(x  )A  is  an  isomorphism  of  Y 

o  o 

means  geometrically  that  our  local  parametrlzatlon  is  valid  for  the  segment  of 
the  solution  path  between  the  points  x  and  x^  closest  to  x^,  where  x_  +  V 
and  x^  +  V  are  tangent  to  the  path.  In  standard  terminology  these  are  the  closest 
limit  points  of  the  path  with  respect  to  the  direction  z^. 


3.  Finite-Dimensional  Approximations 

As  before,  suppose  that  the  Information  (2.1)  is  given  and  that  Rj^(F)  #  i|) . 

Me  turn  now  to  the  formulation  of  suitable  approximate  problems  for  (2.2).  In  many 
applications,  we  have  X  =•  Y  IR;  that  is,  a  particular  component  of  X  is  Identified 
as  a  basic  parameter.  Then  only  the  complementary  component  Y  has  to  be  approxi¬ 
mated.  It  is  useful  to  generalize  this  by  assuming  that  a  splitting 

(3.1)  X  •  Z  ©  T,  dim  T  -  1 
has  been  given  together  with  a  mapping 

(3.2)  Q  €  L(X,Y),  qIZ;  Z  -►  Y  an  isomorphism 


7 


which  relates  Z  with  Y.  Now  a  family  (z,  }  of  finite-dimensional  subspaces  of 

n 

Z,  parametrized  by  some  index  h  >  0,  is  assumed  to  be  given  and  we  set 

(3.3) 

Moreover,  let  Pj^:  Y  -»■  Yj^,  h  >  0,  be  a  family  of  projections  such  that 

(3.4)  lim  llP.y  -  yl  i  =0,  any  y  e  Y. 

h^-0  " 

With  this  our  approximate  problem  is  now  specified  by  the  equations 


(3.5) 


^h^*^  =  ^oh  =  Vo- 


where 


(3.6) 


‘'h^“h=’S.’V 


For  ease  of  discussion,  we  call  the  information  given  in  (3.1)-(3.6)  a  basic  approa 
mat ion  of  our  problem  (2.2). 

As  before,  let  ^^(7^)  “  ^  the  solution  manifold  of  the 

original  problem.  Our  question  is  then  whether  the  solution  manifolds  of  the  dis¬ 
cretizations  (3.5)  approximate  M(y^)  when  h  tends  to  zero.  For  this  analysis, 
we  extend  the  discrete  operators  (3.6)  to  all  of  Wc  X  as  follows: 


(3.7)  Fj^;  W-»  Y,  Fj^(x)  =  (I-Pj^)Qx  +  Pj^(F(x)-y^) ,  x  e  W. 

Clearly,  Fj^  is  of  class  C^(W)  for  h  >  0  and  the  following  properties  can  be 
shown  to  hold: 


(1) 

Fj^(x)  =0  if  and  only  if  x  e  and  “  ^oh’ 

(ii) 

'>\w>S,=  ’'h= 

(3.8) 

(ill) 

ker  DFj^(x)  «=  Xj^; 

(iv) 

DFj^(x)  E  L(X,Y)  is  a  Fredholm  operator  of  index  1; 

(v) 

Pj^DF(x)Xj^  “  Yj^  Implies  that  x  e 
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Any  comparison  of  the  solution  manifolds  of  (2,2)  and  (3.5)  has  to  be  done 
locally.  Hence,  let  e  M(y^)  be  given  and  suppose  that  {V,A,z^}  is  a  local 
parametrization  of  the  manifold  at  x^.  Clearly,  this  parametrization  has  to  relate 
in  a  suitable  way  to  the  basic  discretization  introduced  above.  This  relationship 
may  be  expressed  in  various  forms.  For  the  sake  of  brevity,  we  shall  not  go  into 
details  but  assume  here  simply  the  technical  condition 

(3.9)  1  |Df'jj(x^)Ay|  1  ±  y1  |y|  I  .  y  e  Y,  h  £  (0,h^), 

where  y  >  0  is  independent  of  y  and  h,  and  h^  >  0  is  sufficiently  small.  From 

(3.9)  it  follows  that 

(3.10)  ■  ■'h 

for  h  e  (0,h  )  and  hence,  by  (3.8) (v),  that  x  e  R, (F,). 

o  o  1  h 

The  condition  (3.9)  may  be  enforced  in  many  ways.  For  example,  if  we  are  prepared 
to  restrict  the  discretization  and  parametrization  by  the  choice  V  =  Z,  A  =  A^,  and 
Q  =  A^^Q^,  where  e  L(X,Z)  is  the  natural  projection  corresponding  to  the  splitting 
(3.1),  then  DFj^(x^)A  is  the  identity  on  Y  and  (3.9)  holds  with  y  =  1  for  any  h  >  0. 
On  the  other  hand,  as  noted  earlier,  in  many  applications  we  have  Z  =  Y,  in  which  case 
we  want  to  use  V  =  Z,  A  =  I,  and  Q  =  Q^,  and  (3.9)  turns  out  to  be  a  condition  on  F. 
One  such  case  will  be  discussed  in  the  next  section. 

Let  A^:  Y  -*■  X  be  the  affine  mapping  A^^y  =  x^  +  Ay,  y  e  Y,  and  define 
A^"^^W  c  Y ->■  Y  by  Hj^(y)  =  Fj^(A^y),  y  e  Then,  for  h  e  (0,h^),  the  con¬ 
dition  (3.9)  implies  that  DH^(O)  e  L(Y)  is  an  isomorphism  and  ||DHj^(0)  |  ^  y"^. 

By  the  uniform  boundedness  principle,  the  projections  are  uniformly  bounded  for 

h  >  0.  Using  this  fact,  we  may  show  that  there  exists  a  6  >  0,  Independent  of  h, 
such  that 

(3.11)  l|ltti^(y)  -  DH^(0)||  <  jy  whenever  ||yl'l  <  6. 
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This  permits  the  application  of  the  inverse  function  theorem.  As  noted  in  [s],  the 
standard  proof  of  that  theorem  permits  the  derivation  of  a  Lipschitz  estimate  for 
the  inverse  function.  Using  a  slight  extension  of  the  result  in  [3] ,  we  obtain  here 
the  existence  of  a  unique  C^-function 

(3.12)  Gj^:  |  6y)  <=  Y  ^  B(0,6)  c  Y,  =  (I-Pj^)Qx^, 

such  that 

(3.13)  H^(Gj^(y))  =  y,  Y  e 
and 

Since  6  ,y  are  Independent  of  h  and  i  ^  ^  ^  h  -*■  0,  it  follows  that  there 

exists  some  h.  e  (0»h  )  such  that  0  e  U,  for  h  e  (0,h. ).  Hence,  y.  =  G,  (0) 

J*  o  n  i  n  a 

solves  =  0,  and  from  (3.14)  we  obtain  that  llyj^ll  ^  I  I '^oh^  ^ 

h  e  (0,h^).  Moreover,  we  can  show  that  the  corresponding  points  x^j^  +  x^  +  Ay^  £  X 
satisfy  x^j^  e  Xj^  and  solve  the  approximate  problem  (3.5).  Finally,  it  follows  from 
(3.9)  that  x^j^  E  and  ker  DFj^(x^j^)  0  =  (O). 

With  this  we  have  the  following  result: 

Theorem  3,1:  Let  (2.1)  be  given  and  Rj^(F)  ^  i{>.  Moreover,  suppose  that  a  basic 

approximation  (3.1)-(3.6)  has  been  chosen  and  that  {V,A,z^}  is  a  local  parametrization 

of  M(y^)  at  x^  e  M(y^)  such  that  (3.9)  is  satisfied.  Then,  for  all  sufficiently 

small  h  >  0,  the  approximate  problem  (3.5)  has  a  solution  x  .  such  that  x  e  R,  (F,.) 

on  oh  1  h 

and  lim  x^j^  =  x^.  Furthermore,  for  any  z^j^  e  Xj^,  ^  =  V  fl  Xj^,  it  follows 

h^O 

that  ^ ^h’^oh^  ^  local  parametrization  of  the  solution  manifold 

F^  ^^^yQh^^  ^l^^h^  (3.5)  and  hence  Theorems  2.2  and  2.3  apply. 

This  result  provides  the  existence  of  a  local  segment  of  the  solution  path  of 


the  local  approximate  problem  (3.5)  under  the  same  general  conditions  needed  for 
establishing  the  existence  of  the  solution  curve  of  the  full  problem  (2.2).  However, 
if  error  estimates  are  desired,  then  additional  smoothness  conditions  are  needed  for 
F.  More  specifically,  the  following  extension  of  Theorem  3.1  can  be  proved: 

Theorem  3.2:  Suppose  that  the  conditions  of  Theorem  3.1  hold  and,  in  addition, 

that  DF:  W  ->■  L(X,Y)  is  Lipschitz  continuous  on  bounded  subsets  of  W.  By  Theorem 

2.2,  the  original  problem  (2.2)  has  a  solution  segment  x:  J  c  ]R^  ->■  W  defined  on  an 

open  interval  J  «  0  e  J.  Then  there  exists  a  compact  subinterval  J  c  J, 

o 

0  e  J^,  such  that,  for  sufficiently  small  h  >  0,  the  local  approximate  problem  (3.5) 
has  a  solution  segment  Xj^:  and 

(3.15)  I  lx(t)  -  Xj^(t)|  I  <  C|  I  (I-Pj^)Qx(t)|  I  ,  t  e  J^, 

with  a  constant  C  >  0  which  does  not  depend  on  h  >  0  or  t  e  J^. 

The  proof  uses  a  particular  representation  of  the  two  solution  curves  and  is 
based  on  a  globalized  form  of  the  implicit  function  theorem  given  in  [3], 

4 .  Mildly  Nonlinear  Operators 

As  an  example,  we  apply  our  results  to  the  mildly  nonlinear  operators  considered 
in  [3].  In  line  with  the  comments  in  the  previous  section  about  the  condition  (3.9), 
suppose  that  X  =  Y  x  E  and  denote  by  Q  e  L(X,Y)  the  natural  projection  corresponding 

A  A 

to  this  splitting.  Let  X  be  another  Banach  space,  K  c  L(X,Y)  a  compact  operator, 

A  y 

and  G:  X  ->  X  a  given  nonlinear  mapping  of  class  C  ,  r  ^  1.  With  this  we  assume 
that  our  problem  (2.2)  involves  the  operator  F:  X Y  defined  by  F(x)  =  Qx  +  KG(x) , 

X  e  X.  With  the  notation  x  (y,X),  y  e  Y,  X  e  E,  for  the  vectors  x  e  X,  the 
problem  can  then  be  written  in  the  familiar  form 

(4.1)  y  +  KG(y,X)  “  y^. 


As  in  Section  3,  we  introduce  a  basic  approximation.  Then  it  can  be  shown  that 


for  any  local  parametrization  {V,A,z^}  of  M(y^)  at  e  M(y^),  the  condition 
(3.9)  holds  for  all  sufficiently  small  h  >  0.  Hence,  for  mildly  nonlinear  operators, 
our  basic  assumptions  are  generally  satisfied,  and  therefore  Theorems  3.1  and  3.2 
apply. 

It  may  be  interesting  to  see  what  type  of  local  parametrlzatlons  correspond  to 
the  cases  of  nonsingular  points  and  limit  points  considered  in  the  first  two  parts  of 
[3].  A  point  e  X  is  nonsingular  if  DyF(x^)  =  I  +  KD^G(x^)  is  an 

isomorphism  of  Y  onto  itself.  Hence,  in  this  case,  it  follows  that  for  any 
x^  e  R^(F)  we  have  Y rt ker  DF(x^)  =  (0),  whence  for  a  local  parametrization  we  may 
choose  V  =  Y,  A  the  Identity  on  Y,  and  z^  =  (0,1),  the  X-direction.  This  corres¬ 
ponds  to  the  choice  in  the  first  part  of  [3].  Clearly,  our  results  apply. 

A  point  X  =  (y  ,A  )  e  X  is  a  limit  point  if  D  F(x  )  has  a  one-dimensional 
o  o  o  y  o 

null  space  in  Y  and  D  F(x  )  ^  rge  D  F(x  ).  For  x  e  R,(F),  this  is  equivalent 

AO  y  o  Ox 

with  ker  DF(x  )  c  Y.  Let  u  e  Y  be  such  that  D  F(x  )u  =  0,  I  I u|  I  =  1.  Then 
o  y  o  ♦  •  I  • 

ker  DF(x^)  is  spanned  by  (u,0)  e  X  and  we  have  Y  =  rge  DyF(x^)  ®  span  { u)  ,  and 
D^F(x^)  is  an  isomorphism  of  Y^^  =  rge  DyF(x^)  onto  itself.  It  now  follows  that, 
for  x^  e  Rj^(F),  we  may  introduce  a  local  parametrization  with  V  =  Y^ x  ]r,  ^  (U|0), 

and,  for  Instance,  the  isomorphism  A  e  L(Y,V)  specified  by  Ay  =  (yj^,t^),  where 
y  =  y^^  +  tj^u,  y^  ^  ^  exactly  the  approach  in  the  second  part  of 

[3]  and  again  our  results  apply. 

The  example  discussed  here  is  not  the  only  one  to  which  our  results  apply.  More 
about  this  will  be  presented  elsewhere. 
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